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ABSTRACT 



Context. Hyperluminous infrared galaxies (HLIRG) are the most luminous persistent objects in the Universe. They exhibit extremely 
high star formation rates, and most of them seem to harbour an Active Galactic Nucleus (AGN). They are unique laboratories to 
investigate the most extreme star formation, and its connection to super-massive black hole growth. 

Aims. The AGN and starburst (SB) relative contributions to the total output in these objects is still debated. Our aim is to disentangle 
the AGN and SB emission of a sample of thirteen HLIRG. 

Methods. We have studied the MIR low resolution spectra of a sample of thirteen HLIRG obtained with the Infrared Spectrograph 
on board Spitzer. The 5-8 [im range is an optimal window to detect AGN activity even in a heavily obscured environment. We 
performed a SB/AGN decomposition of the continuum using templates, which has been successfully applied for ULIRG in previous 
works. 

Results. The MIR spectra of all sources is largely dominated by AGN emission. Converting the 6 \im luminosity into IR luminosity, 
we found that ~ 80% of the sample shows an IR output dominated by the AGN emission. However, the SB activity is significant in all 
sources (mean SB contribution ~ 30%), showing star formation rates ~ 300 - 3000 Moyr~' . Using X-ray and MIR data we estimated 
the dust covering factor (CF) of these HLIRG, finding that a significant fraction presents a CF consistent with unity. Along with the 
high X-ray absorption shown by these sources, this suggests that large amounts of dust and gas enshroud the nucleus of these HLIRG, 
as also observed in ULIRG. 

Conclusions. Our results are in agreement with previous studies of the IR SED of HLIRG using radiative transfer models, and we 
find strong evidence that all HLIRG harbour an AGN. Moreover, this work provides further support to the idea that AGN and SB are 
both crucial to understand the properties of HLIRG. Our study of the CF supports the hypothesis that HLIRG can be divided in two 
different populations. 

Key words, galaxies: active - galaxies: starburst - galaxies: evolution - X-rays: galaxies - infrared: galaxies 



1. Introduction 

One of the most important results brought by the Infrared 
Astronomical Satellite (IRAS) was the detection of a new class 
of galaxies where the bulk of their bo lometric emission lies in 
the infrared range (Soif er et al.lll98 4b). This population, named 
luminous infrared galaxies (LIRG), becomes the dominant ex- 
tragalactic population at IR luminosities above 10 u Lq, with a 
space density higher than all other classes of galaxies of co mpa- 
rable bolometric luminosity (see Sanders & Mirabel 1996 for a 
complete review of these objects). At the brightest end of this 
population distribution lie the ultraluminous (ULIRG; Lir > 
10 12 Lq) and the hyperluminous infrared galaxies (HLIRG; 
Lir > 10 13 Lq). 

Ever since the discovery of LIRG two possible origins 
have been suggested for the observed luminosities: Active 
Galactic Nuclei (AGN ) and/or stellar formation in s tarburst (SB) 
episo d es dRieke et alJ Il98(3fc ISoifer et al l Il984at IWright etalJ 
1984 1 ; ISoifer et alJ Il987ak ICarico et al.1 Il990t ICondon et alJ 
19911) . The AGN/SB relative contributions to the IR output 
of this kind of sources shows a clear trend with luminos- 
ity: while LIRG and low-luminosity ULIRG are SB-dominated 
sources, the AGN emission is essential to explain the bolo- 



metric luminosity in high-luminosity ULIRG and HLIRG. The 
fraction of sources harbouring an AGN and its relative con- 
tribution increases with IR luminosity (TVeilleux et alJ 1 1999: 
i Rowan-Robinsonl |200ri iTran et alJ 120011; iFarrah et al.l 12 002a: 
Iveilleux et alJl2002l:lFarrah et aUl200dlNardini et alJl2010h . 

Understanding the ultimate physical mechanism which trig- 
gers those phenomena and, particularly in those sources with 
extreme IR luminosities, their interplay, is crucial to obtain a 
complete view of galaxy evolution. Over the last decade com- 
prehensive observations from X-rays to radio have produced a 
consistent paradigm for local ULIRG. Now is broadly accepted 
that these objects are dusty galaxies where fierce star formation 
processes have been tri ggered by mergers or interactions be- 
tween gas-rich gala xies (Mihos & Hernauis3 ll994l : lGenzel et all 
120011; iDasvra et alj |2006l). Only half of them harbour an AGN, 
which usually is a m inor contributor to the total IR emission (cf. 
lLonsdale et alJl2006l for a complete review). 

ULIRG are rare in the Local Universe (ISoifer et al] fl987b). 
but a large number has been detected instead in deep-IR sur- 
veys, being a fun damental constituent of the high-re dshift 
galaxy populat ion ( Franceschini et al.l 1200 lb lElbaz et al.l [2002; 
Lonsdale etaD 12004 iDaddi et all 120051: fFarrah etaiT 12006: 
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ICaputi etalJl2007t iMagdis et alJl2010L 1201 ll) . It has been pro- 
posed that ULIRG at high redshift coul d be at the ori- 
gin of massive elliptical and SO galaxies dFranceschini et al 



1994; iLillv et alJfl999t Idenzel & Cesarskvl 120001: Idenzel et al 
2001). A large fraction of stars in present-day galaxies 
would have been formed during these evolutionary phases. 
However, the paradigm described above is not so well- 
grounded for the high-z ULIRG population dRodighiero et all 
1201 UlDraper & Ballantvnd2012tlNiemi et al.ll2012l) neither for 
higher luminosity objects like HLIRG. 

HLIRG seem to be composite sources, i.e. AGN and SB 
phen omena are both needed t o fully explain the i r IR emis- 
sion dRowan-Robinsonl 120001: IFarrah et al.1 l2002at iRuiz et alJ 
l2010t iRowan-Robinson & Wandl2010t lHan & Hanll2012h . and 
only about a third has been found in interacting systems 
(Farrah et all l2002bl) . This suggests that most HLIRG are not 
triggered by me rgers, as also found for a significant number of 
high-z ULIRG (ISturm et alJ l2010t iDraoer & Ballantvnel I20H 
lNiemietalJl2012h . If HLIRG and high-z ULIRG share other 
properties, local HLIRG, since they are brighter and easier to 
study with the current astronomical tools, could be key objects 
to understand the high-z population of ULIRG. Moreover, as 
HLIRG could represent the most vigorous stage of galaxy for- 
mation, with SFR > 1000 Mqvt -1 , they are unique laboratories 
to investigate extremely high stellar formation, and its connec- 
tion to super-massive black hole (SMBH) growth. 

The relative contributions of AGN and SB to the bolomet- 
ric luminosity of HLIRG is still in debate. Early studies of 
small samples of HLIRG found contradictory results, with some 
authors sug gesting that the IR emission arises predominantly 
from a SB (Fra ver et alJ 1 199811 19991) with star formation rates 
(SFR) of the order 1000 M0yr _1 , and other au thors suggest- 
ing that HLIRG are powered by a dusty AGN (Granat o et alJ 
1996; lEvans et all 1 19981: lYun & Scoviliel Il998l) . The analysis 
of the IR Spectral Energy Distributions (SED) of larger sam- 
ples of HLIRG r evealed that about half o f these sources are 
AGN-dominated dRowan-Robinsonl 120001: IVerma et al l [2002; 
IRowan-Robinson & Wangll2010l) . Our an alysis of the HLIR G's 
broadband (from radio to X-rays) SED dRuiz et al.l 12010 ) ob- 
tained consistent results, showing that both AGN and SB com- 
ponents are needed to reproduce the full emission of these ob- 
jects, being the bolometric luminosity of the AGN dominant, or 
at least comparable to the SB luminosity, in ~ 70% of the studied 
HLIRG. These studies were however based on samples biased 

to ward AGN. 

IFarrah et alJl2002al (hereafter F02) built a sample of HLIRG 
with a negligible bias toward AGN and it was observed with 
the Submillimetre Common-User Bolometer Array (SCUBA; 
iHolland et alJll999l) . Submillimetre data introduced tight con- 
straints on SB luminosities in the IR SED analysis. They found 
most HLIRG being AGN-dominated sources, but with a signifi- 
cant c ontribution ( > 20% ) due to star formation in all objects. 

In iRuiz eTaD d2007l) we tried to disentangle the AGN and 
SB emission of HLIRG following a different approach. Since 
the 2-10 keV X-ray emissi on above ~ 10 42 e rg s _1 is the "smok- 
ing gun" of AGN activity (Mushotzky 2004), we studied the X- 
ray spectra of a sample of fourteen HLIRG observed by XMM- 
Newton. We found that all the detected sources (ten out of four- 
teen) were AGN-dominated in the X-ray band, but most were 
underluminous in X-rays with respect to the predicted lumi - 
nosity by the standard AGN SED from lRisaliti & Elvis! d2004h . 
given their IR luminosities. We also found strong evidence 
that a significant fraction of these objects is heavily obscured, 
reaching the Compton Thick level (hydrogen column, Nr > 



10 24 cm" 2 ), as suggested by previous studies (Wilmanetal. 
12003b llwasawa et al.ll2005l) . Such absorption level allows, in the 
best cases, only the detection of reflected X-ray emission. To 
avoid the strong obscuration effects and in order to obtain a more 
complete and independent view of these sources, we need to 
move to another spectral window. 

In this paper we present a study of a sample of HLIRG in 
the Mid-Infrared (MIR) band, A = 5 — 8 |om, a spectral range 
very efficient in detecting AGN emission and less affected by 
absorption than X-rays. Several diagnostic methods are avail- 
able to unravel the AGN and SB activity through the MIR spec- 
tra, e.g. the study of high-ionization emission lines and the 
polyc yclic aromati c hydrocarbon (PAH) fea t ures dLaurent et al.l 
2000; lArmus et al.|[2007t ISpoon et al J 120071: IFarrah et al.ll2007l) 
or the analysis of the continuum around 4 [am (Risalit iet al.l 
2006). Blind statistical techniques, like principal component 
analy sis, have been al s o successfully em ployed with MIR spec- 
tra dWang et al.ll201 It iHurlev et afTl2012l) . However, given the 
limited size of our sample, we cannot apply such statistical meth- 
ods. 

The MIR continuum of pure AGN and pure SB show smal l 
dispersion below 10 |jm dNetzer et al.ll2007LlBrandl et al.ll2 006). 
allowing the use of universal templates to reproduce the AGN 
and SB emission in sources where both physical processes are 
present. Hence, we can unravel the AGN and SB components 
of composite sources modelling their MIR spectra with such 
templates. The key reason for using the continuum emission at 
A a 5 - 8 (J.m as a diagnostic is the difference of the 6 |xm- 
to-bolometric ratios between AGN and SB (approximately two 
orders of magnitude larger in the former). This makes the de- 
tection of the AGN component possible even when the AGN is 
heavily obscured and/or bolometrically weak compared to the 
SB. 

SB/AGN continuum spec tral decomposition has been suc- 
cessfully applied i n ULIRG dNardini et al.l 120081 [20091 [20loi: 
iRisahti et al.ll20101) to disentangle the emission of both compo- 
nents. A significant number of H LIRG has been obs erved with 
the Infr ared Spectrograph (IRS, iHouck et all |2004|) on board 
Spitzer dWerner et alj|2004b . obtaining good quality spectra to 
apply this diagnostic technique. 

The outline of the paper is as follows: Section|2]describes the 
selected sample of HLIRG. The IRS data reduction is explained 
in Sect. [3] In Sect. |4] we explain the decomposition process, the 
model adopted and the results obtained. Section [5] discuss the 
results of the MIR spectral analysis and Sect. [6] summarizes our 
conclusions. 

The Wilkinson Microwave Anisotropy Probe ( WMAP) con- 
cordance cosmology has been adopted along the paper: Hq = 70 
km s -1 Mpc 1 , Q m = 0.27, Q A = 0.73 dKomatsu et al.ll2009l) . 

2. The sample 

This paper is based on the analysis of the F02 sample of HLIRG. 
Nine out of ten sources in that sampl43 have been observed with 
the Spitzer-IRS in its low-resolution mode and the data are pub- 
licly available in the Spitzer Archive. 

The F02 sample was originally selected in a manner inde- 
pendent of obscuration, inclination or AGN content. Moreover, 

1 There were eleven sources in the original sample, but we excluded 
one obj ect, IRAS 13279+3401, because the redshift in the literature (z 
= 0.36, IRowan-R obinson 2000) is incorrect, misclassifying the source 
as a HLIRG. We have shown, using recent optical and MIR spectra 
jRuiz et al.ll2010l) . that its redshift is z ~ 0.02. Our estimated IR lumi- 
nosity for this source is hence 3 x 10 10 Lq. 
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its parent sample is composed by objects fou nd from direct op- 
tical f ollow-up of 60 |J.m or 850 |xm surveys dRowan-Robinsonl 
2000) and hence it is statistically homogeneous and complete. 

Given these selection criteria, F02 claim that the sample 
should be entirely free from AGN bias. However, since the par- 
ent sample is composed of sources selected from surveys at 
different wavelenghts, this could produce significant biases in 
favour or a gainst 'hot' dust, which is more comm on in AGN 
than in SB dde Griip et al]|1985b ISoifer et al.lll987al ). Thus, ob- 
jects from 60 |xm surveys are more likely to contain hot dust 
heated by AGN, while those sources from 850 |_im surveys will 
be biased towards cold dust heated by star for mation. However, 
we mu st note that only a small fraction of the iRowan-Robinsonl 
( 120001) sample is selected from 850 [im data (less than 10%), 
and a similar fraction is found in the F02 sample. Moreover, 
the effect of the AGN hot dust in composite sources with equal 
AGN and SB bolometric output is severely diminis hed beyond 
~ 30 rest-frame (cf. Fig. 1 of lNardini et alj2009l) . Only three 
HLIRG in the F02 sample are at z ~ 1, so the effect of hot dust is 
minimal. The potential AGN bias is therefore negligible and we 
consider the F02 sample suitable for drawing global conclusions 
about the HLIRG population. 

Only five sources of this sample have been observed by 
XM M-Newton, an d just two of these five HLIRG were detected 
dRuiz et al.ll2007l) . In order to study the relation between the 
X-ray and MIR emission for this kind of sources, we included 
four additional HLIRG to our sample. These objects were se- 
lected from the X -ray-selected HLIRG sample we analysed in 
iRuiz et all (120071) . They all were detected by XMM-Newton and 
observed with Spitzer-IRS. Three of these additional HLIRG 
were found via compa risons to known AGN (Tables 2-4 of 
IRowan-Robinsonl 12000) and therefore its inclusion leads to a 
non-negligible bias toward AGN content in our final sample. 
Nevertheless, in some ocassions we can derive conclusions about 
the global population of HLIRG restricting our results to the F02 
sample (nine out of ten F02 sources are included in our sample). 

There are thirteen HLIRG in our final sample. Accordingly 
to their optical spectra, six are type 1 AGN (Seyfert 1 or QSO), 
four are type 2 AGN (Seyfert 2 or QS02) and three are SB galax- 
ies (see TableQ]i. The IR SED of all sources have been observed 
and studied using ra diative transfer models (RTM) through sev - 
eral papers (F02; IRowan-Robinsonl 120001: IVerma et al l 120021) . 
AGN and SB components are both needed to reproduce the IR 
emission of these HLIRG, the AGN output being dominant in 
most sources (eight out of thirteen). 

Regarding only to the objects in the F02 sample, five are 
type 1 AGN and four are optically classified as narrow-line (NL) 
objects. All type 1 sources and one NL object show an AGN- 
dominated IR SED, and the remaining (three) NL sources show 
an SB-dominated IR SED, accordingly to the F02 analysis. 

The broad band SED (from radio to X-rays) of most sources 
(nine out of thirteen) have also been studied using observational 
templates for the AGN and SB compo nents, offering co nsistent 
results with those obtained using RTM (IRuiz et al.ll2010h . 

3. Data reduction 

All sources from our HLIRG sample have been observed by the 
IRS onboard Spitzer in the low-resolution mode (see Table [TJ 
and the data were publicly available. All observations were op- 
erated in the standard staring mode using the two low-resolution 
modules (SL and LL). The Basic Calibrated Data (BCD) files, 
bad pixels masks (BMASK) and flux uncertainty files were 
downloaded from the Spitzer Heritage Archive vl .5. All files are 



products of the official Spitzer Science Center pipeline (version 
S18.18). 

The BCD pipeline reduces the raw detector images from in- 
dividual exposures: removes the electronic and optical artifacts 
(e.g. dark current, droop effect, non-linear effects, detection of 
cosmic-ray events), jail-bar pattern and stray light, and perform 
a flat-field correction. BCD files, along with the correspond- 
ing uncertainty images and BMASK files, are reliable inputs for 
SMART, an IDL-ba sed processing and analysis tool for IRS data 
dHigdon et al.ll2004 . 

We processed the BCD pipeline files using SMART 8.2.3. 
We cleaned the BCD files of bad pixels using IRSCLEAN 2.1 
and we combined the individual exposures for a given ExpID 
(uncertainty files were used to perform a weighted-average). In 
order to remove the background and low-level rogue pixels we 
subtracted the two observations in the nodding cycle for each 
module and order (SL1, SL2, LL1, LL2). We checked that all 
sources were point-like and well centered in the slits and finally 
we ran an automatic, tappered column (window adjusted to the 
source extent and scaling with wavelength) extractiorQ to obtain 
the flux and wavelength calibrated spectra. We checked that in- 
dividual orders did match in flux after extraction, so no further 
scaling corrections were needed. 

Two spectra were obtained for each source, one for each nod. 
We calculated an error-weighted average of the two nod spectra^ 
and they were de-redshifted to rest frame. Figure Q] shows the 
final spectra. 

4. SB/AGN spectral decomposition 

4. 1 . Model for the MIR emission of HLIRG 

Thanks to the unprecedent sensitiviy of Spitzer-IRS, MIR stud- 
ies of large samples of luminous AGN and SB have found a sig- 
nificant spectral homogene i ty in the SEP of the two separate 
classes (Brandl et all l2006t INetzer et all I2007I: iMullanev et alj 
I201 lb iHer nan-Caballero & Hatziminaoglou; 201 ll) . For high lu- 
minosity sources, this is particularly true within the 5-8 [im 
range, showin g nearly constants shapes of both the PAH com- 
plex in SB dBrandl et all I2006I) and the continuum in AGN 
(INetzer et alJl2007t) . 

This small spectral dispersion allows the use of universal 
AGN/SB templates to parameterize the observed energy output 
of objects where both phenomena are present, like ULIRG or 
HLIRG, in the wavelength range below 8 |xm. This technique 
has been extensively tested in ULIRG, showing that the observed 
differences of the PAH equivalent width and strength are caused, 
in most sources, by the r elative contribution of the AGN contin- 
uum and its obscuration (iNardini et al.ll2008l I2009L 120 1 Oh . 

The bolometric luminosities of HLIRG and ULIRG differ 
no more than about an order of magnitude (even less for the SB 

2 We also teste d the advanced optimal (PSF weighted) method 
iLebouteiller et al.ll2010h to extract the spectra. Since all sources are 
bright and point-like, the results from both extraction methods were 
similar. However, we found that the spectra obtained through the ad- 
vanced optimal method show a strong sinusoidal artifact in most SL2 
spectra. We com pared our spectra with th e Cornell Atlas of Spitzer/IRS 
Sources (CASIS, Lebouteiller etaU201lD and we found no similar fea- 
tures at those wavelengths. We chose then the tapered column extraction 
to perform the subsequent analysis. 

3 Averaging the spectra introduces a systematic error related to the 
flux differences between the two nod spectra. This difference is usually 
due to the presence of other sources in the slit affecting the background 
subtraction, or to pointing errors result ing in a slight shift in the disper- 
sion direction ( Lebouteiller et al. 201 1). 
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component). Moreover, a significant fraction o f HLIRG seems t o 
share common properties with local ULIRG dRuiz et al . 2010). 
Thus, we can rely with r eason able confidence on the model pro- 
posed by Nar dini et al.l d2008l hereafter N08) to reproduce the 
5-8 ii.m emission in HLIRG. The observed flux at wavelength A, 
f° hs , is given by: 

f obs _ Ant r ,AGN„-T(/t) , / 1 \ „SBl n x 

where a(, is the AGN contribution to (the intrinsic flux den- 
sity at 6 |am) while m^ gn and are the AGN and SB emission 
normalized at 6 [im. This model has only three free parame- 
ters: a^, T(, (the optical depth at 6 urn) and (the flux normal- 
ization). Additional high-ionization emission lines and molecu- 
lar absorption features (due to ices and aliphatic hydrocarbons), 
whenever present, were fitted using Gaussian profiles. 

Starburst emission 

As stated above, starbursts of different luminosities present a 
small spectral dispersion in their MIR emission below 8 |xm. 
Given the large number of physical variables (e.g. the initial 
mass function, the duration and evolution of the burst, dust grain 
properties) involved in determining the observational proper- 
ties of an SB, it i s in fact surprising finding such similarity. 
Bran dl et~aT] (|2006) suggested that this homogeneity can derive 
from the spatial integration of unresolved star-forming spots. We 
can then assume that an observational template will accurately 
reproduce the SB emission of our sources in the 5-8 \xm range. 
This emission is characterized by two prominent PAH features 
at 6.2 and 7.7 [im (Fig.|2]i. 

To increase the robustness of our analysis we have em- 
ployed four different SB templates from N08 (a) and from 
iHernan-Caballero & Hatziminaogloull201 ll (b, c, d): a) the av- 
erage MIR spectrum from the five brightest ULIRG among pure 
SB sources in the N08 sample; b) the average spectrum of ob- 
jects optically classified as SB or HII in NED and with vL v (7 
|j.m)< 10 44 ergs -1 (16 sources); c) the average spectrum of SB- 
dominated objects in the MIR (257 sources) and d) the average 
spectrum of objects classified as ULIRG (184 sources). As ex- 
pected, all SB templates show similar feature^ and small dis- 
persion below 8 ii.m (see Fig. [2j. 

AGN emission 

In our range of interest the AGN spectrum is dominated by the 
cooling of small dust grains that are transiently heated up to tem- 
peratures close to the sublimation limit. This process is expected 
to produce a nearly featureless power law cont inuum. However, 
PAH features are found in s everal quasa r s (|Schweitzer et alj 
2006; Lutz et alJl2007l l2008h . iNetzer et alJ (l2007h adscribe the 
FIR emission of a set of luminous PG QSO to cold dust in ex- 
tended regions, thus probably of stellar origin. After removing 
this galactic contamination, the average spectrum can be mod- 
eled with a single power law from ~ 3 [im to the 9.7 \im silicate 
bump. Since the spectral dispersion within this range is small 
for pure luminous QSO, we can model the AGN emission in our 
sources with the same power law with a fixed spectral index: 
fx °c /l 8 . Again, to increase the robustness of our results, we 

4 Slight differences in the shape of the lines and their intensity ra- 
tio are probably caused by variations of the obscuration in the origi- 
nal sources: higher lumi nosity sources usually show higher obscuration 
dRigopoulou et alll999h . 



alternatively fitted all sources allowing the slope of the power 
law to vary in the range 0.2-1.8. This interval covers the disper- 
sion found in the spe ctral index for the AGN emission in ULIRG 
dNardini et alJl2010h . 

Our model includes an exponential attenuation in the AGN 
emission to take into account the reddening of the NIR radia- 
tion due to any compact absorber in the line of sig ht. The optica l 
depth follows the conventional law t(A) oc A' 1 - 75 (lDraind ll989). 
The SB emission does not need a similar correction because the 
possibles effects of obscuration are already accounted in the ob- 
servational templates. 

4.2. Results 

The model presented in the pre vious section was implemented 
in Sherpa (IFreeman et alJl200ll) . the modelling and fitting tool 
includ ed in the software package CIAO 4.3 dFruscione et al.l 
2006). The 5-8 spectra were fitted using a Levenberg- 
Marquardt y 2 minimization algorithm. The errors of the pa- 
rameters (within lcr confidence limit) were calculated through 
Montecarlo simulations. 

One object was barely detected by Spitzer and in two addi- 
tional sources our model was not suitable (see Sect. 14.31 for a 
further discussion on these objects), but the spectra of the re- 
maining ten sources were well reproduced by our model (see 
Fig.[3]andTable|2l). 

We found no significant difference between the best-fit mod- 
els using the four SB templates. All templates offer a similar 
X 2 and consistent values for the model parameters. We have 
therefore included in Table [2] only the results obtained using 
the N08 SB template. However, we should note that the use 
of the ULIRG template provides systematically lower values 
of <X(, in all sources. This result is probably caused by the 
presence of an AGN component in the template, as expected 
in an average ULIRG spectrum including high-luminosity ob- 
jects (Nardi ni etalJ l2010: H ernan-Caballero & H atzimina ogloul 
I20T1 . 

Allowing the spectral index to vary did not enhance the fit in 
most sources. A statistically significant improvement (in terms 
of F-test) was found only in three sources with high optical depth 
(beside IRAS 18216+6418, which is discussed below). However 
we must note that our model is severely degenerated in the op- 
tical depth-spectral index plane (see Fig. 0). The shape of the 
5-8 [im continuum cannot be well established in sources show- 
ing high absorption. Hence, if the spectral index is a free pa- 
rameter, the minimization algorithm decreases the optical depth 
while increasing the spectral index. We selected the model with 
fixed spectral index as our best fit in those sources showing ob- 
scuration features outside the 5-8 |xm range (e.g. 9.7 \im silicate 
absorption). 

The 6 |J.m emission is largely dominated by the AGN com- 
ponent (ff6 > 0.9), which also dominates over the whole 5-8 |xm 
spectral range. IRAS F00235+1024 and IRAS F16124+3241 are 
the only two sources where the SB and AGN components are 
comparable, but the AGN emission is still above the SB. The 
PAH emission at ~ 6.2 and ~ 7.7 \im in these two sources is well 
reproduced by our SB template, suggesting that our assumption, 
i.e. the MIR SB emission of HLIRG can be modelled by the MIR 
emission of SB-dominated ULIRG, is correct. 

All sources optically classified as type I AGN show mi- 
nor absorption (except IRAS F10026+4949), while all sources 
classified as NL or type II AGN show high absorption, (except 
IRAS 12514+1027, which shows a strong 9.7 [im silicate ab- 
sorption feature). 
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In two sources (IRAS 09104+4109 and IRAS 12514+1027) 
we detected a spectral feature at ~ 7.7 ii.m that can be modelled 
as a gaussian emission line. In both cases the line is statistically 
significant (in terms of F-test), although it is not needed to obtain 
a good fit. The line can be interpreted as the unresolved [NeVI] 
fine-structure line at 7.65 \xm. Due to its high ionization energy 
(126 eV), it seems to be connected to the hardest nuclear activity. 
This emission featu re has been detected in previous s tudies of 
local ULIRG dNardini et alJl2009t IVeiUeux et alJ2009h . 

4.3. Notes on particular sources 
IRAS 00182-7112 

This Compton-Thick QS02 shows a heavily absorbed MIR 
spectrum. At first sight, our model seems to reproduce well 
the spectrum above 6.5 iim (see Fig. |3(a)| i if two absorption 
lines at 6.8 and 7.3 \xm are included (C -H stretching mod e 
of hydrogenated amorphous carbon (HAC: lFurton et alJll999lV ). 
However, the lack of the 6.2 (j.m PAH emission feature and the 
overall shape of the whole MIR spectrum (see Fig. |1 (a)[ > strongly 
suggest that th e 7.7 \xm feature is actually not associated with 
PAH emission. iSpoon et alJ d2004a) suggest that the 5.7-7.7 |xm 
range is dominated by a broad absorption complex due to water 
ice and hydrocarbons. Thus, our simple model cannot reproduce 
the complex spectra of this object (as suggested by the poor fit 
we obtain, with a reduced^- 2 ~ 20). 

The spectrum shows no signs of SB activity within the 5- 
8 i_im range, since the PAH features are totally suppressed, and it 
is very similar to the MIR spectra of deeply obscured AGN (e.g. 
NGC 4418, NGC 1068). This result, along additional evidence 
from other wavelengths (e.g. strong iron K„ emis sion line with ~ 
1 keV equivalent width in the X-ray spectrum. lRuiz et aU l2007). 
suggests that the power source hiding behind the optically thick 
mate rial could be a buri ed AGN. 

ISpoon et"aT1 (|2004a) present a detailed study of the whole 
IRS spectrum: the strong absorption and weak emission features 
in the 4-27 \±m spectrum suggest the existence of a dense warm 
gas cloud close to the nucleus of the source. Based on the strengh 
of the 11.2 \im PAH feature, there is also evidence of star for- 
mation activity away from the absorbing region, that could be 
responsible for up to 30% of the IR luminosity of the system. 
For further analysis and comparison we have employed this es- 
timate of the SB relative contribution. 

IRAS 09104+4109 

This is a radio loud QS02 and the MIR spectrum is probably 
dominated by synchrotron radiation, modifying the shape of the 
AGN emission. Although our model offers an excellent fit in a 
pure statistical sense (y 2 /d.o.f. 0.6), the selected AGN tem- 
plate just models the reprocessed emission by dust. Hence, in 
this case, the physics predicted by our model could be inaccu- 
rate. 

IRAS 18216+6418 

The continuum of this QSO seems to be flatter than the adopted 
power law. This is the only non-obscured source where a free 
spectral index model offers a clear improvement in the best-fit. 
The flatter continuum could be just the natural dispersion around 
the template: a flatter AGN s lope has been detected in ~ 10% of 
ULIRG (Nardini et ai] |2010l) . However, we must note that this 
object, although classified as a radio-quiet AGN, show many 



properties of radio-loud sources (F02). Thus we cannot reject 
the posibility that the flatter continuum is caused by synchrotron 
contamination. 

IRAS F23569-0341 

This source was barely detected by Spitzer and the extracted 
spectrum had a very poor S/N ratio, so it was rejected in the 
subsequent analysis. This object has not been detected neither 
with ISO or SCUBA (F02). 

5. Discussion 

The 5-8 \xm spectra of most ULIRG show signatures of AGN ac- 
tivity, but the relative contribution of the AGN emission at 6 (xm 
spans a broad range, fro m complete SB-dom i nated to com plete 
AGN-dominated output dNardini et al.ll2008l I20091 l2010h . Our 
spectral decomposition, on the other hand, clearly states that 5- 
8 \xm spectra of HLIRG are strongly dominated by the nuclear 
AGN emission reprocessed by the dusty torus. This is consis- 
tent with previous studies showing that the fraction of ULIRG 
harbouring an AGN and the relativ e contribution of this com- 
ponent increase wi th IR luminosity (IVeilleux et al.l fl999. 2002; 
iNardinietallkolOl) . 

This result is also in agreement with previous stud- 
ies in other wavelen gths which found AGN activity in 
most HLIRG (F02: IVermaet alJ l2002t iRowan-Robinsonl 
2000; iRowan-Robinson & Wang|l2010t iRuiz et al.ll2007l I2010). 
Moreover, based on the completeness and non AGN bias of the 
F02's sample, our analysis is a strong direct evidence that all 
HLIRG harbour an AGN. 

5.1. Star formation rate from PAH emission 

Many studies have shown the tight correlation between the 
presence of Polycyclic Aromatic Hydrocarbon (PAH) features 
in the MIR spectrum and the presence of starburst activity 
dGenzel et al.lll998l:lRigopoulou et alJl2000HBrandl et~aiT2 006). 
The PAH emission arises in the photo-dissociation region that 
lies between the HII region of an SB and the surrounding molec- 
ular cloud where the stars are formed. The PAH spectral features 
are, in addition, very uniform, especially below 8 um (cf. Fig. 1 
of N08 or Fig. 2 from lSargsvan & Weedmanll2009l) . 

The 5-8 |xm spectra of our HLIRG are strongly dominated 
by the AGN emission, so we cannot obtain a reliable direct mea- 
sure of the PAH emission. Instead, we estimated the emission of 
this spectral feature through the SB component obtained with 
our spectral decomposition. For each source, we calculated a 
rough estimate of the PAH peak flux at 7.7 [am using the best-fit 
model0 

/f (7.7 urn) = /** (1 - a 6 ) uf (7.7 jim). (2) 

After converting these fluxes into luminosities, we estimated 
the star formation rate (SFR ) using the relation obtained by 
ISargsvan & Weedmanl(l2009l) : 

log SFRpah = log [ALf (7.7 urn)] - 42.57 ± 0.2, (3) 

where SFRpah is the star formation rate in solar masses per year 
and /IL^ B (7.7 (j.m) is the PAH luminosity at 7.7 \±m in ergs per s. 

5 w;j B (7.7 (i.m) was estimated interpolating the N08 SB template at 
7.7 p.m. 
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Our estimate of the 7.7 \im flux includes both the PAH emis- 
sion and the underlying SB continuum due to dust emission, but 
Eq. [3] is calibrated respect to the total 7.7 |xm flux, including 
also both components. In any cas e, the continuum contributio n 
is typically just ~ 10% of the total (Sargsyan & Weedman 2009). 
Since is close to unity in all our sources, the relative error of 
the SB contribution at 6 [im is high. Thus our SFR estimates 
cannot be accurate beyond the order of magnitude. 

Figure [5] compares our SFR estimates wi th those obtained 
through modelling the IR SED using RTM dRowan-Robinsonl 
2000; IVerma et alj2002l F02). Both techniques seem to find con- 
sistent results. Despite most being AGN-dominated objects (see 
Sect. [572] below), we still found high star forming activity, with 
SFR ~ 300 - 3000 M© yr" 1 (see Table©. 

Figure [6] shows the SFR we calculated versus (a) total and 
(b) SB IR luminosity. We estimated the SB luminosity as: 



Lfi =(l-a IR )L n 



(4) 



where am is the fractional contribution of the AGN to the 
total IR luminosity (see Sect. 15. 2\ and Lm is the total IR 
emission (8-1000 jam) c alculated through the IRAS fluxes 
(ISanders & Mirabel 1 19961) . When only an upper limit is avail- 
able for the IRAS fluxes at 12 or 25 |xm, we have simply assumed 
the one half of this limit. This convention provides a fair aprox- 
imation of the actual flux densities even for sources at moderate 
redshift dNardini et alJl2^09ll20Toh . 

We performed several statistical tests to check any possible 
correlation between the SFR and the IR luminosities. Since our 
estimates of and SFR are both dependent of the parameter 
(see Eqs. |2]andj2l), we used a generalize d Kendall's Tau partial 
correlation tes£j(Akritas & Siebert 1996). The result shows that 
both quantities are correlated with a probability ~ 99.99% > 3<x. 
This tight correlation between the star formation and SB IR lu- 
minosity is expected from sev eral theoretical results ( Ke nnicuttl 
ll998tlDraine & Lill200lll2007L and references therein). 

In the case of total IR luminosity versus SFR we can ap- 
ply a generalized Kendall's Tau tesjj, since the estimate of the 
IR luminosity is completely independent of our estimate of the 
SFR. They show a weaker, but slightly significant, correlation 
(the probability is 98% > 2<x). This slight correlation and the 
high SFR suggest that SB emission could be a significant con- 
tributor to the IR output (see Sect. 15.21 ). 



5.2. AGN contribution to the IR luminosity 
We define the 6 ^m-to-IR bolometric ratio as: 



R 



(5) 



where Fir is the total IR flux (8-1000 |xm) estimated as in 
ISanders & Mirabell (Il996l) . The hypothesis that the IR and to- 
tal bolometric fluxes are almost coincident is typically a fair as- 
sumption in ULIRG, but in principle we cannot adopt this hy- 
pothesis in our sources. The study of the HLIRG's broad band 
SED has revealed that, for an important number of sources clas- 
sified as HLIRG, a non-negligible fraction of the bolometric 



6 Kendall's Tau test is a non-parametric correlatio n test based on the 
Kendall rank coefficient (Kendall & Gibb ons! 1 1990h . The generalized 
version is an extension to include censored data (e.g. upper limits). We 
used the code developed by the Center for Astrostatistics to perform this 
test (http : //astrostatistics . psu . edu/statcodes/cens_tau I. 

7 We used the ASURV software for this test dlsobe et aljl986h . 



luminosity is emitted outside the IR range dRuiz et al.1 12010). 
Therefore we limited our analysis to the IR luminosity. 

We can derive a connection between R and the parameter org: 



R = 



^AGN ^SB 



a 6 R SB + (1 -a 6 ) R AGN ' 



(6) 



where R ACN and R SB are the equivalents of R for pure (unob- 
scured) AGN and pure SB, as defined in Eq.|5] 

R and ac are both known quantities. Assuming a^ as the in- 
dependent variable and R as the dependent variable, we can fit 
Eq. [6] to the values observed in our sources, considering R ACN 
andfl SB as free parameters. We found R AGN = 0.36 +0.08, R SB = 
(8 + 3) x 10 -3 . In spite of the limited size of our sample and 
the narrow range for q-6, these results are in good agreement 
with those obtained by N08 for a larger sample of ULIRG: 
^agn _ o.32+^, J R SB = (ll.7^)xl0 -3 . Fig.[7]shows/? versus 
Qf6 for our sample of HLIRG, including our best-fit and the result 
from N08. 

Using a 6 , R AGN and R SB (using the N08 values for the last 
two quantities), we could estimate the fractional AGN contribu- 
tion (am) to the IR output of each source: 



7AGN 



ClR 



CC6 



Fir a 6 + CR AGN /fl SB )(l-a 6 )' 



(7) 



Figure [8] compares <?ir estimated through our MIR spectral 
analysis and through the analysis of the IR SED using RTM 
dRowan-Robinsonl feoOOt IVerma et all 120021 F02;). Both esti- 
mates seem to be consistent in most sources. Letting aside those 
sources where our model is not suitable (IRAS 09104+4109 and 
IRAS 00182-71 12), we found that the IR luminosity of seven out 
of ten HLIRG is AGN dominated. However, the SB contribution 
is also significant for all sources, spanning from ~ 20% - 70% 
(see Table 0. Using only the sources from the F02 sample, the 
mean SB contribution is ~ 30%, close to 35% obtained by F02. 
This significant SB contribution is consistent with the correla- 
tion we found between SFR and IR luminosity and the high SFR 
estimated (see Sect. 15. U . Our analysis confirms the idea that star 
formation and accretion into SMBH are both crucial phenomena 
to explain the properties of these extreme objects. 

5.3. Covering Factor 

An important physical parameter to unveil the distribution of 
dust in the nuclear environment of an AGN is the dust covering 
factor (CF), i.e. the fraction of sky covered by dust viewed from 
the central engine of the AGN. The CF is also critical to under- 
stand the fraction of direct nuclear emission that is absorbed by 
dust and re-emitted in the IR range and hence has a significant 
effect in the bolometric luminosity corrections of AGN. 

The ratio between the thermal AGN luminosity (i.e., the IR 
reprocessed emission due to heated dust), ^j™, and the primary 
AGN luminosity (above ~ 1 [am, i.e., the accretion disk bolomet- 
ric emission), Lp™ , is commonly interpreted as an estimate of 
the dust C F (Maioli no et ai1l2007[ iRowan-Robinson et al.ll2009b 
Roseboom et al. 2012|): 



CF(dust) 



j AGN 

TH 
j AGN ' 

BOL 



(8) 



The thermal AGN luminosity can be estimated through the 
MIR continuum emission. Assuming that the IR AGN luminos- 
ity is dominated by the dust emission and using the best-fit pa- 
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rameters obtained in our MIR spectral decomposition, the ther- 
mal emission is computed as: 



j AGN 



^AGN 



(9) 



X-ray emission is a primary product of the central engine 
of AGN, so X-ray luminosity is usually a good proxy to es- 
timate the bolometric AGN luminosity (Maiol ino et al.1 120071 
IRowan-Robinson et al.l2009l). We appl ied the bolometric correc- 
tion estimated by Marconi et alJ ([2004): 



log 



j AGN 
BOL 



= 1.54 + 0.24X + 0.012X 2 - 0.0015X 3 



(10) 



where X = log^Lg™/^©) _ 12 and Lx is the intrinsic (i.e., 
absorption-corrected) 2-10 keV luminosity. 

Therefore we can estimate the CF of those sources for which 
MIR and X-ray data are both available. Eight out of thirteen 
HLIRG in our sample ful fill this condition . The X-ray luminosi- 
ties were obtained from iRuiz etail 120071 IRAS F00235+1024 
and IRA S 07380-2342 are pro bably Compton-Thick (CT) 
sources (IRuiz et al.l 120071 l2010h so the upper limits on their 
X-ray luminosities have been increased by a factor of 60, the 
average ratio betw een intrinsic and o bserved X-ray luminosi- 
ties in CT sources dPanessa et al.1 12006). IRAS 14026+4341 was 
not detected in our X- ray analysis, but it has a counterpart in 
the 2XMMi catalogue dWatson et al.ll2009l) . We estimated its 2- 
10 keV luminosity using the 2XMMi X-ray fluxes. Since thi s 
object seems to be an X-ray absorbed QSO dRuiz et al.l 12010). 
we applied the correction for CT objects described above to cal- 
culate its intrinsic X-ray luminosity. 

Table [4] shows all the derived AGN luminosities and final 
estimates of the dust CF. The covering factors versus the AGN 
luminosities at 6 [J.m are plotted in Fig. [9] For comparison, the 
plot includes the average CF estimated for local QSO0 

Four out of eight objects show a dust CF consistent with ~ 1 
(at 2<x level) and sistematically above the average CF of local 
QSO (~ 0.5). This result and the heavy X-ray absorption shown 
by most of these sources (all but one source, IRAS 14026+4341, 
show signatures of CT absorption in their X-ray emission, see 
Wilman et al. 2003; Iwasawa et al. 20051 IRuiz et alJ2007ll20Tol 



Vign ali et alj|201 lb point towards large amounts of gas and dust 



enshrouding t heir nuclear environment as it has been found in 
local ULIRG dSpoon et al.ll2004bl: IVerma et alJl200l lYan et alJ 
l2010h . 

We found one source with CF ~ 0.2. This source, 
IRAS 18216+6418, is a QSO with no sign of X-ray or MIR 
obscuration. The low CF is consistent with previous studies 
finding a decrease in the CF of QSO with in creasing lumi- 
nosity dMaiolino et al.ll2007l: iTreister et ai1l2008l). The decrease 
could be explained by "receding toru s" models dLawrencdl99U 
iMaiolino et al]|2007l: lHasingeril2008l) : low-luminosity AGN are 
surrounded by a dust torus of obscuring material covering a 
large fraction of the central source. High-luminosity AGN would 
be able to clean out the environment ionizing the surrounding 
medium or blowing it away through outflowing winds. Hence 
the opening angle of the torus would be larger and the covered 
solid angle would be lower (assuming the height of the torus is 
not luminosity-dependent). 



8 The average CF w as e stimated through a di rect integration of the 
iRichards et ajT lT2006) and iHopkins et alT j2007h average AGN SED. 
The direct AGN emission was calculated integrating the SED at wave- 
lengths below 1 (im, and the thermal emission integrating the SED at 
wavelengths above 2 \im. 



The remaining three sources show CF >> 1. A dust CF 
abov e unity is usually interpreted as evidence of CT obscura- 
tion (Ro wan-Robinson et al.l 120091) . As stated above, two ob- 
jects (IRAS F00 235+1024 and IR AS 07380-2342) were not de- 
tected in X-rays (Ruiz et al.l2007l) and their broadband SED sug- 
gest that their low X-ray emission is due to high obscuration 

T II —j 

(Ruiz et al. 2010). Hence the bolometric AGN emission is prob- 
ably underestimated in these sources. Observational and theoret- 
ical studies of CT AGN predict even larger values (i.e., >> 60) 
for the intrinsic-to-observed X-ray luminosity ratio, depending 
on the amount of absorption and on the viewing angle with re- 
spect to the obscuring torus dHaardt et al.l ll994; Iwasawa et al. 
1997). The dust CF of these sources would be consistent with 
unity if the largest correction factor (~ 1000) were applied. 

The third remaining source, IRAS F12509+3 122, is a QSO 
showing an X-ray emission significantly lower than the one 
predicted for a high-luminosity QSO, given its IR luminosity 
(Rui z et al.ll2.Q10l) . This suggests that its bolometric AGN emis- 
sion is also underestimated. However, the low X-ray luminosity 
of thi s object cannot be related to X-ray absorption (Ruiz et al. 
2007). Alternatively, an increase with Eddington ratio has been 
observed in the X-ray bolometric correction of AGN (Aem - 
L B OL/L E dd), from k x ~ 15 - 30 for X Edi £ 0.1 to K X ~ 70 - 120 
for A Edd > 0.2 (IVasudevan & F abian 2009]). Us ing the SPSS op- 
tical sp ectrum of IRAS F12509+3 122 and the iMcLure & Jarvisl 
(2002) relation between black hole masses and Mgll emission 
line widths, we estimated A Edd ~ 0.5 for this object. This high 
Eddington ratio suggests that the X-ray bolometric correction for 
this source could be higher than that obtained in Eq. [10] Further 
studies are needed to investigate how the Eddington ratio influ- 
ences the X-ray luminosity of HLIRG. 

6. Conclusions 

We studied low resolution MIR spectra of thirteen HLIRG ob- 
served with Spitzer, nine of them also observed with XMM- 
Newton. We modelled their 5-8 \xm spectra, using the AGN/SB 
spectral decomposition technique developed by N08, and we es- 
timated the contribution of each component to the total IR lumi- 
nosity. 

We found that all HLIRG in the sample harbour an AGN that 
clearly dominates the 5-8 \m\ spectrum. Given the completeness 
of the F02 sample (nine out of ten F02 objects are included in 
our sample), this is strong evidence suggesting that all HLIRG 
harbour an AGN. We have also found that the AGN component 
seems to be the dominant power source of the total IR output 
in most sources. However, all sources, even the AGN-dominated 
HLIRG, show a significant SB activity, with a mean SB contri- 
bution of ~ 30%. The SFR, estimated through the PAH emission, 
is also very high in all sources (~ 300 - 3000 Mgyr -1 ). 

These results are in agreement with previous stud- 
ies of HLIRG using theoretical models to reproduc e their 
IR SED (F02: iRowan-Robinsonl 120001: IVerma et all 120021 
IRowan-Robinson & Wangl 1201 Ol) . providing further support to 
the assumptions of these detailed models. The AGN relative con- 
tribution to the total IR output is also consistent with those es- 
timated through SED modeling. The large mean AGN contribu- 
tion we found (~ 70%) is consistent with previous studies of 
ULIRG, pointing towa rds an increase of the AGN emission wit h 
increasing luminosity ( Veilleux et al. 2002; Nar dini et al.l l2010). 
Our work confirms the crucial role of both AGN and SB phe- 
nomena to explain the properties of these extreme sources. 

Using X-ray and MIR data we were able to estimate the cov- 
ering factor (CF) of these HLIRG, finding that a significant frac- 
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tion (seven out of eight) have CF> 1. Most of these sources 
with large CF also show heavy X-ray absorption and high op- 
tical depth or absorption features in their MIR spectrum. This 
strongly suggests that the nuclear environment of these sources 
is heavily enshrouded by large amounts of gas and dust, as it has 
been observed in local ULIRG. 

F02 proposed that HLIRG could be divided in two pop- 
ulations: (1) mergers between gas-rich galaxies, as found in 
the ULIRG population, and (2) young active galaxies going 
through the ir maximal star fo rmation periods whilst harbouring 
an AGN. In lRuiz et alJ (l2010) we fo und strong ev i denc e support- 
ing the two-population hypothesis. lHan & Hani d2012l) studied 
the same HLIRG sample using AGN and SB theoretical models 
and bayesian techniques to fit their IR SED. They found further 
evidence that the physical properties of AGN and SB harboured 
by HLIRG could be divided into two populations: SB-dominated 
HLIRG (corresponding to our "ULIRG-like" objects) show a 
higher fraction of OB stars and the starburst region is more 
compact, while AGN-dominated objects show a dustier tho- 
rns. Additiona l evidence favouring this tw o-populations idea is 
presented by iDraper & Ballantvnd d2012l) . They suggest, based 
on recent Herschel observations, that major mergers cannot ex- 
plain the total population of luminous high-z ULIRG harbouring 
AGN. 

All sources showing CF ~ 1 (except I RAS 14026+4341) 
were classified as "ULIRG-like" objects in iRuiz et all (2010) 
while IRAS 18216+6418, with a significantly lower CF, 
was classified as a "non-ULIRG" source. Furthermore, 
IRASF00235+1024 and IRAS F15307+3252, both sources 
showing a high CF, were observed by the Hubble Space 
Telescope (HST) and both present signs of recent mergers. 
On the other hand, IRAS 18216+6418, al so observed by HST , 
shows no signs of interactions or mergers dFarrah et alj|2002bl) . 
The study of the dust CF adds another piece of evidence favour- 
ing that HLIRG can be separated into two differentiated popula- 
tions, although further studies with larger samples of HLIRG are 
needed to obtain stronger conclusions. 
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Table 1. Spitzer observations of HLIRG. 



Sources " 


Sample b 


Type c 


R.A. 


DEC. 


z 


AOR" 


Exposure (s) 


Date e 




IRAS 00182-7112 


X 


QS02 


00 20 34.7 


-70 55 27 


0.327 


7556352 


94.4 


2003-11 


14 


IRAS F00235+1024f 


F,X 


NL-SB 


00 26 06.5 


+ 10 41 32 


0.575 


12237056 


94.4 


2005-08 


11 


IRAS 07380-2342f 


F, X 


NL 


07 40 09.8 


-23 49 58 


0.292 


12236032 


12.6 


2005-03 


23 


IRAS 09104+4109 


X 


QS02 


09 13 45.4 


+40 56 28 


0.442 


6619136 


73.4 


2003-11 


29 


IRAS F 10026+4949 


F 


Syl 


10 05 52.9 


+49 34 42 


1.120 


4738560 


14.7 


2004-04 


17 


IRAS F 12509+3 122 


FX 


QSO 


12 53 17.6 


+31 05 50 


0.780 


12236800 


62.9 


2005-06 


06 


IRAS 12514+1027 


X 


Sy2 


12 54 00.8 


+ 10 11 12 


0.32 


4978432 


94.4 


2005-02 


06 


IRAS 14026+4341 1 


F,X 


QSO 


14 04 38.8 


+43 27 07 


0.323 


4374016 


102.8 


2005-05 


22 


IRAS F15307+3252 


X 


QS02 


15 32 44.0 


+32 42 47 


0.926 


4983552 


243.8 


2004-03 


04 


IRAS F16124+3241 


F 


NL 


16 14 22.1 


+32 34 04 


0.71 


4984832 


243.8 


2005-03 


19 


ELAIS J 1640+41 


F 


QSO 


16 40 10.1 


+41 05 22 


1.099 


11345664 


243.8 


2005-08 


13 


IRAS 18216+6418 


F, X 


QSO 


18 21 57.3 


+64 20 36 


0.297 


4676096 


31.5 


2004-04 


17 


IRAS F23569-0341 


F 


NL 


23 59 33.6 


-03 25 13 


0.59 


12235776 


243.8 


2004-12 


13 



Notes. (a) Sources marked with f have been o bserved with XMM-Newton, but not detected. (6) F: sources from the F02 
sample; X: sources from the lRuizetafl d2007h X-ray-selected sample. (c> Optical spectrum classification. (< " Astronomical 
Observation Request (AOR) of the Spitzer observation. <c) Date of the Spitzer observation. 



Table 2. Best-fit model parameters (c.f. Sect.0). 



Sources " 


Model 6 


x z /d.oI. 


/■"' (mly) c 




Spectral index e 


t 6 > 


Additional features 


IRAS F00235+1024 


A 


238/75 


86 + 5 


0.982 ± 0.002 


0.8 


2.57 ± 0.06 






B 


218/74 


46+f 


969+ 012 


> 1.0 


1 q^+o.53 
^ -0.05 




IRAS 07380-2342 


A 


22/76 


213 ± 12 


0.988 + 0.006 


0.8 


0.49 ± 0.05 






B 


16/75 


130 + f 


0.990 ± 0.009 


1 64 +0M 


<0.3 




IRAS 09104+4109* 


A 


44/69 


176 ±3 


0.995 ± 0.002 


0.8 


0.69 ± 0.02 


[NeVI] emission line g 




B 


26/68 


104 ±11 


0.996 ± 0.004 


1.65 ±0.17 


0.16 ±0.10 


[NeVI] emission line^ 


IRAS F 10026+4949 


A 


19/88 




> 0.988 


0.8 


63 +0 06 

-0J2 

6 +(M 






B 


19/87 


150 ± 60 


> 0.992 


7 +0 - 7 

" -0.5 




IRASF12509+3122 


A 


51/85 


25.3+S;* 


QSQ+0.005 

• -0.003 


0.8 


< 0.04 






B 


51/84 


25 ±6 


962+™ 

-0.003 


0.87!^ 


< 0.2 




IRAS 12514+1027 


A 


116/73 


57.4 + 1.0 


0.958 ± 0.003 


0.8 


< 0.02 


[NeVI] emission line- e 




B 


114/72 


57 ±7 


0.956 ± 0.004 


77+ 01 

O. / 1 _ 012 


< 0.12 


[NeVI] emission line 8 


IRAS 14026+4341 


A 


33/75 


75 ±2 


0.977 + 0.004 


0.8 


0.06 ± 0.02 






B 


31/74 


70 + | 3 


0.978 ± 0.004 


93 +0 02 

-0.31 


< 0.18 




IRAS F15307+3252 


A 


238/84 


38 + 2 


0.981 + 0.004 


0.8 


0.64 ± 0.04 






B 


173 / 83 


21^ 


0.975 ± 0.006 


> 1.4 


07 +0 24 

"•"'-0.03 




IRAS F16124+3241 


A 


38/78 


6.6 ± 1.3 


0.8735? 


0.8 


1.3 ±0.2 






B 


38/77 


7±4 


0-8^ 


1 3 +0 - 5 




ELAIS J 1640+41 


A 


33/88 


17.3 + 1.1 


0.993 + 0.006 


0.8 


0.13 ± 0.06 






B 


32/87 


15 +6 . 

136 r 22 

159' 


> 0.988 


1 06 +002 

1,uu -0.60 


< 0.4 




IRAS 18216+6418 


A 


117/77 


> 0.999 


0.8 


< 0.001 






A 


18/76 


n OOI +0 005 

v - yyl -n ma 


<0.4 


< 0.17 





Notes. <n) Sources where our model is not reliable are marked with *. { - b) HLIRG emission model (A: AGN power law slope 
fixed at 0.8; B: AGN power law slope allowed to vary between 0.2 and 1.8). See Sect. l4.1l for a complete discussion of the adopted 
model. The best-fit model is shown in boldface. (r) Intrinsic flux at 6 (im (corrected by dust absorption). Relative contribution 
of the AGN component at 6 \im. {e) Spectral index of the AGN power law. Optical depth at 6 \im of the obscuring material. 
(s) Emission line at 7.65 \im. See Sect l4.2l for a discussion about these spectral features. 
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Table 3. SFR and AGN relative contributions to the IR output, and IR luminosities (in cgs units). 



Sources " 


SFR" 




j — 71 — 

MR 


j fJB<? 

[xlO 46 ] 




[M yr'] 




[xlO 46 ] 


IRAS F00235+1024 


1 500 

uuu 1(x)0 


0.68 + 0.12 


20 ±8 


6 + 4 


IRAS 07380-2342 


< 700 


°- 82 S 


14.2 + 0.6 


3 + 2 


IRAS 09104+4109* 


< 600 


Q9+0.08 
-0.07 


13.8 ± 1.3 


1 ' 1 -09 


1KAS rlUUzo+494y 


< DUU 


> U.5 


1 AC\ i ACt 

14U ± 41) 


< 10 


IRAS F12509+3 122 


2000 + 1400 


049+o.i4 
53 + ° : ^ 

"■-"-0.11 


36 ± 13 


18 + 8 


IRAS 12514+1027 


400 ± 300 


3.7 + 0.4 


l. 7 +0.7 


IRAS 14026+4341 


400 + 300 


60+ 012 
-Q!Q 


6.2 + 0.5 


2 5 +0 ^ 

-0,7 


IRAS F15307+3252 


3000 + 2000 


59 +!i " 


52+ 17 


21 +i6 

~ 9 « 


IRAS F16124+3241 


i3 00+ i7oo 
< 1200 


0.25!; 


16 + 6 


12+ 6 


ELAIS J 1640+41 


0.82!;« 


< 120 


< 20 


IRAS 18216+6418 


300+300 


0.80+n,- 


8.0+1.4 


1 6 +L0 



Notes. (n) Sources where our model is not reliable are marked with *. {b) Star formation rates estimated through 
PAH intensity (see Sect l5.lt . {c> AGN contribution to the total IR luminosity (see Sect. 15.2b . {d> Total IR luminosity 
(8-1000 um), estimated through IRAS fluxes. (c) IR luminosity (8-1000 urn) of the SB component (see Sect.[5~Tb. 



Table 4. AGN luminosities (in cgs units) and CF values. 



Sources 


L x " 


/ 1MR b 

^AGN 


A L A (6 //m) c 


j TH d 

^AGN 


CF f 




[xlO 44 ] 


[XlO 46 ] 


[xlO 46 ] 


[xlO 46 ] 




IRAS F00235+1024f 


< 1.6 


< 0.6 


6.1 ± 1.2 




> 30 


IRAS 07380-2342f 


< 2.0 


<0.8 


3.0 + 0.3 


n 


> 11 


IRAS 09104+4109 


20!f 




6.8 + 0.3 






IRASF12509+3122 


1.8 + 0.2 


0.73 ± 0.09 


3.9 + 0.3 


12! 4 


16 +6 


IRAS 12514+1027 


2 +4 5 


04 +L41 

"•"^-0.03 


0.88 ± 0.06 


T 7+1.0 
^•'-0.7 


7^+3000 


IRAS 14026+434 It 


3.0+1.8 


1 4+1.0 

-0,9 


1.27 + 0.07 


4 +L4 
-1.0. 


2.8^ 


IRAS F15307+3252 


31+ 8 


34 +9 


8.1+0.7 


25 + » 


0.7 + 0.3 


IRAS 18216+6418 


37 + 4 


43 + 5 


2.11+0.10 


6.6^ 


15 +006 

u - 1J -0 04 



Notes. (n) Absorption-corrected X-ray luminosity, from iRuiz et alJ ([2007). A factor of 60 iPanessa et al.l I2006T) 
has been applied to transform the observed X-ray luminosity to intrinsic X-ray luminosity to the sources marked 
with f . <6> Bolometric AGN luminosity (i.e., intrinsic AGN emission above ~ 1 p.m). {c) Absorption-corrected AGN 
luminosity at 6 p.m. id) Reprocessed AGN luminosity (i.e., thermal emission due to heated dust). (f,) Covering factor. 
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Wavelength (/xra) 




5 10 15 

Wavelength (/xra) 




Wavelength (/xra) 




Wavelength (/xm) Wavelength (/i-m) Wavelength (/xm) 




Fig. 1. MIR rest-frame spectra of HLIRG obtained with the Infrared Spectrograph onboard Spitzer. The grey shaded area is the lcr 
uncertainty region. The vertical dashed lines limit the spectral region where our decomposition technique was applied. 
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6 8 10 12 14 6 8 10 12 14 

Wavelength (|xm) Wavelength (|j,m) 




6 8 10 12 14 6 8 10 12 14 



Wavelength (|J,m) Wavelength (|J,m) 

Fig. 2. Average SB spectra used as temp lates in our model: Template (a) is from N08; templates (b), (c) and (d) from 
iHernan-Caballero & Hatziminaogloul d201 ll) (see Sect. |4]for details). The grey shaded area is the estimated Icr dispersion. 
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IRAS00182(x /d.o.f. : 19.6; SB: 1) 



IRASF00235 (///d.o.f. : 3.2; SB: 1) 



IRAS07380 (x"/ d.o.f. : 0.3; SB: 1) 
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Wavelength (u m) 
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Wavelength (\im) 




5.5 6 6.5 7 

Wavelength (u m) 



IRAS09104 (/ /d.o.f. : 0.6; SB: 1) 



IRASF10026 (///d.o.f. : 0.2; SB: 1) 



IRASF 12509 (/// d.o.f. : 0.6; SB: 1) 




6 6.5 7 

Wavelength (u m) 




5.5 6 6.5 7 7.5 

Wavelength (u m) 




5.5 6 6.5 7 7.5 

Wavelength (u m) 



IRAS12514(/ /d.o.f. : 1.6; SB: 1) 



IRAS14026(x / d.o.f. : 0.4; SB: 1) 



IRASF15307 (x I d.o.f. : 2.8; SB: 1) 




5 5.5 6 6.5 7 7.5 

Wavelength (u m) 




5.5 6 6.5 7 7.5 

Wavelength (um) 



0.008 
' OiKK. 
ODO 




5.5 6 6.5 7 

Wavelength (um) 



IRASF16124 (x 2 / d.o.f. : 0.5; SB: 1) ELAISJ164010 (/ 2 / d.o.f. : 0.4; SB: 1) IRAS18216 (j I d.o.f. : 0.2; SB: 1) 




5 5.5 6 6.5 7 7.5 8 5 5.5 6 6.5 7 7.5 8 5 5.5 6 6.5 7 7.5 8 



Wavelength (u m) Wavelength (u m) Wavelength (u m) 



Fig. 3. 5-8 |J.m spectra and best fit models (blue solid line). The AGN component is ploted as a red dashed line. For clarity, the SB 
component (green dotted line) is ploted only on those sources where the AGN and SB contributions are comparable. 
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Fig. 4. Spectral Index-Optical depth confidence regions for IRAS 07380-2342 (left) and IRAS F16124+3241 (right). The contours 
represent lcr, 2<x and 3<x confidence levels. 
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1000 
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Fig. 5 . Comparison of the SFR estimated through PAH emission and through IR SED modelling using RTM (Far rahetalJ l2002a: 
lRowan-Robinsodl2000tlVerma et al.ll2002l) . 
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Fig. 6. SFR estimated using the PAH emission at 7.7 |im versus (a) total and (b) starburst IR luminosities. The open circles mark 
those sources where our model is inaccurate. 
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Fig. 8. Comparison of the AGN contribution to the IR luminosity estimated through MIR spectral decomposition and through RTM. 
The blue squares are type I AGN and the red triangles are type II AGN and SB. The open circles mark those sources where our 
model is inaccurate. MIR contribution for IRAS 00182-71 12 from lSpoon et al.1 d2004al) . 
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Fig. 9. CF versus the AGN luminosity at 6 \xm. The blue dashed line is the estimated average CF for local QSO (see Sect. 15.3b . 
Symbols as in Fig. [8] 
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